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ABSTRACT 

We present photometric redshifts and spectral energy distribution (SED) classifications for a sample 
of 1542 optically identified sources detected with XMM in the COSMOS field. Our template fitting 
classifies 46 sources as stars and 464 as non-active galaxies, while the remaining 1032 require templates 
with an AGN contribution. High accuracy in the derived photometric redshifts was accomplished as 
the result of 1) photometry in up to 30 bands with high significance detections, 2) a new set of SED 
templates including 18 hybrids covering the far-UV to mid-infrared, which have been constructed by 
the combination of AGN and non-active galaxies templates, and 3) multi-epoch observations that have 
been used to correct for variability (most important for type 1 AGN) . The reliability of the photometric 
redshifts is evaluated using the sub-sample of 442 sources with measured spectroscopic redshifts. We 
achieved an accuracy of (JAz/{i+z,p^^) = 0.014 for <22.5 (a'A2/(i-i-ZspEe) ~ 0.015 for i^^ <24.5). 
The high accuracies were accomplished for both type 2 (where the SED is often dominated by the host 
galaxy) and type 1 AGN and QSOs out to z = 4.5. The number of outliers is a large improvement 
over previous photometric redshift estimates for X-ray selected sources (4.0% and 4.8% outliers for 
'^*AB <22.5 and i^^ <24.5, respectively). We show that the intermediate band photometry is vital to 
achieving accurate photometric redshifts for AGN, whereas the broad SED coverage provided by mid 
infrared {Spitzer/IUAC) bands is important to reduce the number of outliers for normal galaxies. 
Subject headings: AGN: general, photometric redshift 
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1. INTRODUCTION 

It is now well established that AGN play a signifi- 
cant role in cosmic evolution - tracing the growth of 
super massive black holes in galaxy nuclei and pro- 
viding important energetic feedback to the host galax- 
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ies and the IGM. Characterization of the AGN prop- 
erties and their redshifts is therefore of paramount im- 
portance to understanding cosmological evolution. The 
mass of the central super massive black hole (SMBH) 
directly correlates with the mass and t he velocity dis- 
persi o n of the bulge of its host galaxy (jCebhardt et al.1 
l2000t iFerrarese fc MerritH l2000f) and these correlations 
apparently pe rsist both for normal, le ss-active galaxies 
and for AGN (jMcLure fc Dunlop|l2003 ). This is strongly 
suggesting a co-evolution of galaxies and their nuclear 
BHs. Unfortunately, given the low fraction of galaxies 
with high luminosity AGN and the fact that those with 
high luminosity may give a biased view of the overall 
importance of the AGN feedback, the co-dependent evo- 
lution of galaxies and BH has been pursued very little 
at high redshift, where the majority of both galaxy and 
AGN evolution occurs. Here we address this problem - 
attempting to identify a much larger sample of AGN over 
a range of luminosity and to provide classifications and 
redshifts for much fainter samples. 

At present, no criterion exists to identify all AGN in 
a field using a single photometric band. Optical, in- 
frared, radio, narrow bands and X-ray techniques each 
yield incomplete, partially overlapping, sub-samples of 
the overall AGN population. A complete census can 
only be achieved by combining the results of the var- 
ious selection techniques (Zamorani et al., in prepara- 
tion). Nevertheless, the X-ray band appears to be very 
efficient and rather comp lete for most classes of AGN 
(|Brandt fc Hasingeil[200a ). 

The COSMOS survey provides full multi-wavelength 
data-sets from X-ray to radio a nd we focus here on 
XMM-COSMOS selected sources (|Hasinger et al.ll200 iD 
using spectral energy distributions (SEDs) to disentan- 
gle between non-active galaxies and AGN and to provide 
accurate photometric redshifts. Photometric redshifts 
for normal (inactive) galaxies can reach an accuracy 
of crAz/(i+z) = 0.05 ( Grazian et al.l l2006t lllbert etldl 
[20061) or better (|Gabasch et all [2004: Wolf et al.l[200l 
and most recently better than 0.02 (Ilbert et al. 2008). 
In contrast, fo r AGN the most a ccurate results have 
been 0.06-0.1 (I Zheng et al.l l2004t iPolletta et all [20071 : 
iRowan-Robinson et al.ll200S ) for redshift z<l. 

Uncertainties in the SED fitting of AGN come from 
two major sources. First, AGN can vary with time and 
quasi-simultaneous photometry is required to obtain a 
snapshot SED. Second, the observed emission is often a 
superposition of AGN and host galaxy. Depending on 
the properties of the host and the relative power of the 
nucleus, the different spectral bands can be dominated 
by either of the two components. For example, NGC3079 
is a Seyfert 2 galaxy, identified on the basis of its optical 
nucleus and a radio jet, yet its Spitzer/IRS spectrum is 
often used as a te mplate for non active, s tarburst galax- 
ies in the infrared (jWeedman et al.ll2006af ). On the other 
hand, Mrk231, a Seyfert 1 and ultra-luminous infrared 
galaxy (ULIRG), exhibits strong silicate absorption in 
the IRS spectrum and is therefore highly obscured with 
much AGN (and starburst) luminosity re-emitted in the 
far infrared. Clearly, to analyze samples which, as XMM- 
COSMOS, contain both galaxies and AGN, one needs: 1) 
a fine sampling and full spectral coverage with photome- 
try spanning from the UV to mid-infrared and 2) spectral 



energy distribution templates representing the full range 
of AGN and QSOs together with variable mixing frac- 
tions of AGN and host galaxy emission. 

Only recently, with Spitzer and GALEX providing full 
SED templates of A GN beyond the optical bands (e.g. 
IPolletta et al.l [2007[ l such analysis has become possible 
at moderately high redshift. 

In this paper we derive photometric redshifts for the 
AMM-COSMOS sources - mostly dominated by an AGN 
(see § Ej) - achieving an accuracy comparable that ob- 
tained for non activ e galaxies. The COSMOS field 
(jScoville et al.l l2007bl ) with its size (2 square degrees) 
and deep photometry in 30^ bands from the radio to 
the X-ray, is an excellent testbed for such a census. 

In § [21 we describe the sample and the photometric 
data set. The variability analysis is presented §[3| In §|3] 
the procedure used for the AGN classification and pho- 
tometric redshift determination is discussed. The results 
are given in § [S] and the description on the released cata- 
log are described in § [5] The discussion and the conclu- 
sions end the paper in §[7] and §[5| Throughout this work 
we use AB magnitudes and assume Ho=70kms"iMpc~\ 
r2A=0.7, and VLm^Q-^- 

2. THE DATA SET 

2.1. The XMM-COSMOS sample 

The COSMOS field has been o bserved with XMM- 
Newton for a total of 1.5 5 Ms ([Hasinger et al.l l2007l : 
ICap pelluti et all [20071 12008D at the homogeneous depth 
of ~ 50 ks. The final catalogue includes point-like 
sources detected above a given threshold with a max- 
imum likelihood detection algorithm in at least one of 
the soft (0.5-2 keV), hard (2-10 keV) or ultra-hard (5-10 
keV) bands down to limiting fluxes of 5x10"-^^, 3xl0~^^ 
and 5xl0~^^ erg cm~^ s~^, respectively (see Cappel- 
luti et al. 2007, 2008 for more details). The adopted 
threshold {Lik > 10) corresponds to a probability of 
~ 4.5 X 10~^ that a source is a background fluctua- 
tion. A detailed X-ray to optical association using COS- 
MOS optical, near- and mid-infrared catalogs was per- 
formed using the Maxim um Likelihood (ML) technique 
(|Brusa et al.l[2007l[200l . For the sub-field of the XMM- 
COSMOS covered also by Chandra, the optical identi- 
fications have been augmented with the more accurate 
Chandra positions (Civano et al. 2008; Elvis et al. 2008). 
As a result, the sample can be divided into the following 
subgroups: 

• 1528 sources (78.9%) have unique and secure opti- 
cal counterparts within 5."0 from the AMM coordi- 
nates. We will refer to these sources as the "secure 
sample" . 

• 52 sources are brighter than i* = 16.5 and are sat- 
urated in many (if not all) optical bands. 37 of 
these are associated with bright stars, and 1 is a 
galaxy with spectroscopic redshift at 0.122. These 
38 sources are removed from the input catalog as 
all the photometry is saturated. We will refer to 
the remaining 14 sources as the "bright sample". 

• 10 sources are undetected in the optical catalog 
and the identification of the counterpart is done in 
the near- and mid-infrared using the K-band (Mc- 
Cracken et al. 2008, in preparation) and IRAC 
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catalogs^^. For this "opticaUy undetected sample" 
not more than 5 photometric points are available. 

• 233 sources have an ambiguous optical cross- 
identification. For these objects, the ML technique 
yielded 2 sources with similar probability. This 
sample is referred to as having "ambiguous coun- 
terparts" . 

The photometric redshift determination is tuned using 
only the "secure" and the "bright" samples (1528+14 
= 1542 sources); all of the results and the discussion be- 
low refer only to them. The optically undetected sample 
will be discussed separately in § 15. 3| while for the sample 
with ambiguous counterparts, we will provide photomet- 
ric redshift for the primary and secondary counterpart 
in the catalog but they are not included in the redshift 
distribution or in the discussion. 

2.2. Multi-wavelength data 

For the 1542 XMM sources associated with an opti- 
cal/infrared counterpart, the photometric redshift has 
been computed using broad and intermediate band pho- 
tometry in the range 1350A< A <8/im using GAL EX, 
Subaru, SDSS, CFHT and Spitzer/IKKG. 

The optical/near -infrared photometry is extracted 
from the recently updated COSMOS photometric cat- 
alog (Capak et al. 2008) . Compared to the previous 
catalog (Cap ak^t al.l 120071 ). the new photometry imple- 
ments: 1) a more accurate algorithm for source detection, 
2) new, deeper (and with better seeing) u*- and Kg-band 
imaging from CFHT and 3) new J-band data (UKIRT) 
and 12 additional intermediate and 2 narrow bands (from 
Subaru; Taniguchi et al 2007, Taniguchi et al. 2008, in 
preparation). The fluxes were measured in fixed aper- 
tures of 3."0 diameter, on PSF-matched i mages (FWHM 
of l."5). Simulations (jCapak et al.ll20d7l ) show that the 
fiux in 3."0 aperture corresponds to 0.759 of the total 
flux for a point-like source. 

For each optical/infrared counterpart of a source in 
the XMM-COSMOS sample, the closest GALEX source 
from the de-blended, P SF-fitted GALEX-CO'$MO^ cata- 
log of lZamoiski et al.l (|2007t l was considered for inclusion 
in the photometry data set. Only 625 sources do have a 
GALEX counterpart. 

IRAC photometry was also extracted from the IRAC- 
COSMOS catalog. In this case, the possible IRAC coun- 
terpart of each XMM source was taken to be the clos- 
est IRAC source within l."0 of the optical counterpart. 
For approximately 5% of the XMM-COSMOS sources, 
blending of multiple IRAC sources makes it impossible 
to reliably extract a correct IRAC flux. Additional 5% 
fall in masked regions in the IRAC images and therefore 
lack mid-infrared photometry. 

The 3."8 aperture fluxes given in the COSMOS IRAC 
catalog (aperture 2 in the catalog) where then trans- 
lated to total fluxes by division by factors of 0.765, 0.740, 
0.625, 0.580 for 3.6^m, 4.5^m, 5.8/im, and 8/im, respec- 
tively (see the readme file associated to the catalog and 
Surace et al. 2004 for details). 

Lastly, we multiply the total fiux from GALEX and 
IRAC by 0.759, in order to re-scale the photometry to the 

http://irsa.ipac.caltcch.edu/data/COSMOS/ 



optical and near-infrared one. All photometry has been 
transformed to the AB magnitude system. The available 
bands, depth and observing epochs are summarized in 
Table [J 

3. VARIABILITY ANALYSIS 

AGN exhibit non-periodic fl ux variations on time 
scales of minutes to decades (see lWold et al.|[2007l for a 
review). Thus, photometry collected non-simultaneously 
will be affected by variability and may be not represen- 
tative of a snapshot SED. Wolf et al. (2004) pointed out 
that ignoring variability can significantly increase the er- 
rors in photometric redshift estimations of AGN. 

The COSMOS optical photometry has been acquired 
in five epochs distributed over ~6 years (see column 7 in 
Table [T]). Within each epoch (apart from 2005), several 
individual filter observations were distributed over less 
than 3 months, covering the whole optical range. The 
photometry in each filter has generally a very high cali- 
bration accuracy. This allows to study and correct AGN 
time variability over time scales of years, while shorter 
variability can not easily be addressed. 

Comparing the extracted fiuxes at different epochs, a 
number of objects show clear time variability. As an ex- 
ample, in Figure [1] we show the variability of source XID 
17 and the correction that we apply. In the upper panel 
the raw optical photometry of this source is shown as 
black data points with different symbols for each epoch. 
There is obviously a large scatter between the photom- 
etry at different epochs. The time variability correction 
proceeded in the following way: 




aoo2 



2004 
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2006 



Fig. 1. — Time variability for the source with XID 17. Bottom: 
The source was observed in 5 epochs, each epoch marked with a 
different symbol. The 5 groups are, from left to right: SDSS, Sub- 
aru broad band images (Subarugs)i CFHT i band, the first set 
of intermediate bands Subaru images (Subaru/si) and the second 
epoch of intermediate bands Subaru images (Subaru/B2)- The op- 
tical photometry varied up to 0.6 mag over the 6 years of COSMOS 
data acquisition. The red line connects the median values for each 
group of observations. Top: the photometric points of the source 
before (black) and after (green) our correction. The red and blue 
lines represent the smoothed spectral shape of the object before 
and after the correction for variability (see text for more details). 
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TABLE 1 

Photometric coverage and depth of the COSMOS data-set used for the photometric redshift estimation. 



Filter 


telescope 


effective A 

A 


FWHM 

A 




Depth 
mag^s 


Time of observations 
UTC 


u* 


CFHT 


3911.0 


538 


0.05 


26.50 


Tan 2003 - Anr 2007 


R T 


^1 1 ri!l I'll 


4439.6 


807 


-0.24 


27.00 


Tan 9004 


V 7 


O Li UcLL LI 


5448.9 


935 


-0.09 


26.60 


Feb 2004 


ff+ 

6 


^1 1 rinvn 

O Li L^OjJ. L1 


4728.3 


1163 


0.02 


27.00 


Tan 2005 


r+ 


O LlUdl LI. 


6231.8 


1349 


0.00 


26.80 


Tan 9004 


i+ 


O LllJcll LI. 


7629.1 


1489 


0.02 


26.20 


Tan 9004 

(J (Xll ^UVyrt 


z+ 


^1 1 V^iivn 

O Li L^OjI. LL 


9021.6 


955 


-0.04 


25.20 


Tan 9004 

U (Xll ^\J\J^ 


J 


UKIRT 


12444.1 


1558 


0.12 


23.70 


Mar 900fi 


Ko 


CHFT 


21 480 2 


3250 


-0.05 


23.70 


Mar 2007 

i.VX(ILl f 




CHFT 


7629 


1460 


-0.01 


24.00 


Tan 9004 

(J (Xll .^VJL/t: 




SDSS 


3623 


557 


0.00 


22.00 


Till 2001 

U Lil -L^ W L/ i 


S 


SDSS 


4677 


1171 


0.00 


22.20 


Jul 2001 




SDSS 


6159 


1123 


0.00 


22.00 


Till 2001 

U Li.1 jLj V J 17 i 




SDSS 


7480 


1291 


0.00 


21.30 


Tnl 9001 

■J Lii .^VJL/1- 




SDSS 


8884 


943 


0.00 


20.50 


Till 2001 

U Lil i*WL/l- 


IA427 


^1 1 r^iivn 
O LL L^OjI. ll 


4256.3 


207 


0.04 


25.82 


Tan 900fi 

U (Xll 


IA464 


O Li UcLL Ll 


4633.3 


218 


0.01 


25.65 


Feb 2006 


IA484 


O LlUdl Ll. 


4845.9 


229 


0.00 


25.60 


Tan 9007 

•J (Xll jiUL/ 1 


IA505 


^1 1 riiivn 

O LL L/CiL Li. 


5060.7 


231 


0.00 


25.55 


Feb 2006 


IA527 


^1 1 riiivn 

O LL L^Oji. Li. 


5258.9 


242 


0.03 


25.62 


Tan 2007 

U (Xll -L^ W C 


IA574 


^1 1 Vinvn 

O l-i IjCJLL Ll 


5762.1 


272 


0.08 


25.61 


Tan 2007 

(J (Xll ^VJL/ ( 


IA624 


O LllJcll Li. 




301 


0.00 


25.60 


Dpr 900R 


IA679 


Subaru 


6778.8 


336 


-0.18 


25.60 


Feb 2006 


IA709 


Subaru 


7070.7 


316 


-0.02 


25.65 


Feb 2007 


IA738 


Subaru 


7358.7 


324 


0.02 


25.60 


Jan 2007 


IA767 


Subaru 


7681.2 


364 


0.04 


25.60 


Mar 2007 


IA827 


Subaru 


8240.9 


344 


-0.02 


25.39 


Jan 2006 


IRACl 


Spitzer 


35262.5 


7412 


0.00 


23.90 


Jan 2006 


IRAC2 


Spitzer 


44606.7 


10113 


0.00 


23.30 


Jan 2006 


IRAC3 


Spitzer 


56764.4 


13499 


0.01 


21.30 


Jan 2006 


IRAC4 


Spitzer 


77030.1 


28397 


-0.17 


21.00 


Jan 2006 


FUV 


GALEX 


1551.3 


231 


0.31 


25.69 


Feb 2004 


NUV 


GALEX 


2306.5 


789 


-0.02 


25.99 


Feb 2004 



In order to flatten out the overall spectral shape 
to first order, the wavelength dependence of the 
optical magnitudes was simply smoothed by a run- 
ning average filter. The result of this smoothing is 
shown as the red solid line in the upper panel of 
Figure [TJ 

This smoothed spectral shape was then subtracted 
from each individual data point in order to calcu- 
late the relative magnitude differences shown in the 
lower panel of Figure [1] Here the magnitude differ- 
ences Am are plotted against the observation date 
of the 5 separate epochs, n, of contemporaneous ob- 
servations. For a constant source all points in this 
panel should be within the measurement uncertain- 
ties around Am=0. In this particular example, the 
source XID 17 apparently brightened by more than 
0.6 magnitude over the 6 years. 

For the epochs n— 1,2,4 and 5 the median Am„ 
of the magnitude differences of all optical bands 
observed in this epoch was calculated. For epoch 
n=3, which consists of only one data point, Ama 
was computed by a linear interpolation between 
n=2 and n=4. 

The correction offset for a given epoch n is the dif- 
ference offsetn — Am„ — Ato„=4, where n=4 is 



the epoch of the Subaru/si and IRAC observations 
(Dec 2005-Mar 2006). Each point in the plot rep- 
resents the deviation of the photometry in a given 
filter observed in epoch n, from m„. 

• In order to derive the closest approximation to a 
snapshot SED, we corrected for variability sub- 
tracting the value of offsetn from each of the 
multi-epoch photometry data points, resulting in 
the green data points in the upper panel of Fig- 
ure [1] 

• Similar to the red curve, the thin blue solid line in 
the upper panel of Figure [1] shows the smoothed 
spectral shape of the variability corrected magni- 
tudes. 

In the case of source XID 17 the variability correction 
clearly resulted in a significant reduction of the scatter in 
the photometry and thus an improvement in the photo- 
metric redshift determination. As discussed above, this 
correction scheme only works for slow time variability 
and there are examples of sources, where a significant 
time variability could not be reduced by this method. 
However, given the temporal sampling of the diverse data 
sets, it is the best that could be done and this correction 
was found to significantly reduce the number of outliers 
for the spectroscopic sub-sample (see § ^ . 
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As anticipated from the known variabihty time scales 
of AGN, the AMM-COSMOS sources exhibited a variety 
of behaviors - some sources evolved gradually while oth- 
ers displayed more extreme outbursts. We quantify the 
observed variability with the parameter, 



VAR = 




(1) 



In Figure [2] the distribution of VAR for all the sources 
(black solid line) is shown. 



100 




VAR 



Fig. 2. — Histogram of the parameter VAR for all the sources 
(black solid line) extended (red) and point-like sources (blue). The 
vertical dotted line at VAR=f).2b defines our adopted separation 
between host (low VAR) and AGN (high VAR) dominated sources. 



4. FITTING TECHNIQUE 

To derive the intrinsic SEDs and redshifts of the ob- 
served sources, we used the publicly available LePhare 
code (Arnouts, Ilbert)^^. For this fitting, a library of ex- 
pected intrinsic rest-frame SEDs of the source population 
is supplied and a minimization is used to solve for the 
most likely redshift, SED type and intrinsic extinction. 
A user-supplied prior can also be used. In the following 
we describe the transformation of the template SEDs to 
the observed photometric bands, the selection of tem- 
plate SEDs, the adopted extinction law and the control 
of source variability which is important for AGN. 

The accuracy of the photometric redshift fitting was 
evaluated using a sub-sample of 396 XMM sources with 
i* <22.5 and secure spectroscopic redshift. 90 sources 
had optical sp ectrosco py from Sloan Digital Sky Sur- 
vey ()York et al . 2000; Adelman-McC arth y et ah ,2008.) 
while 306 were observed within the COSM OS collabora- 
tion (235 observed with IM ACS and MMT (jTrump et all 
[2007t iPrescott et all [20061): 7 1 observed by zCOSMOS 
with VIMOS (i Lillv et al.ll200 /)). The spectroscopic sam- 
ple was used to train/evaluate the fitting procedure 

^ = 1.48 X 



median(|zphot - ^spcc|/(l + ^spoc) (|Hoaglin et al.l ll983f) . 
where Zphot and Zspec are the photometric and spectro- 
scopic redshifts, respectively. As the reliability of pho- 
tometric redshifts trained on local/bright spectroscopic 
samples may not hold for low- luminosity or high redshift 
sources for which the SED templates may not longer be 
representative, we have also performed an a posteriori 
test using spectroscopy which became available only re- 
cently for 46 optically faint, K-band and Spitzer/yUFS 
24/im selected sources (Lilly et al. 2008; Kartaltepe et al. 
2008). These additional spectroscopic redshifts indicated 
a high success rate with correctly recovered redshift also 
for these faint sources (see 



4.1. 



Transformation of Template SEDs to Observed 
Photometric System 



and the accuracy quantified as crAz/(i- 



http:/ /www. oamp.fr/ people/arnouts/LE_PHARE.html 



To reduce the spectra template to the same photomet- 
ric system of the observed sources, all SEDs were con- 
volved with the filter transmission curves, which also in- 
clude the instrument efhciency. The same zero-points de- 
rived for the new COSMOS photometric redshift catalog 
of non X-ray detected galaxies (Ilbert et al. 2008) were 
used for the AMM sources. This zero-point calibration is 
extensively explained in Ilbert et al. (2006) and reported 
in column 5 of Tab. [TJ The convolved template SEDs 
were then redshifted in Az= 0.02 wide bins from z=0.02 
to z=:8 and parabolically interpolated for higher redshift 
accuracy. The opacity of the inter-galactic mediu m was 
taken into account as described by Mada^ ()1995( ). 

4.2. Galaxy and AGN Templates 

Ideally, a template library should represent all possi- 
ble SEDs. However, the larger the number of templates, 
the higher the risk of degeneracy, which then can lead to 
an increasing number of catastrophic failures. To iden- 
tify a suitable library, we run LePhare on the spectro- 
scopic sample of 396 sources with i^^ <22.5. Here, we 
fixed the redshifts to the spectroscopic values and let the 
code select the best fitting SED template. Initially, the 
entire set of 28 templates from SWIRE^^, was used, in- 
cluding elliptical and spiral galaxies, starbursts, ULIRG, 
low luminosity AGN and QSOs (QSOl, QS02, TQSOl, 
BQSOl etc; see PoUetta et al. 2007 for details on these 
templates). We also added two composite low- and high- 
luminosity QSO templates from SDSS and one template 
of a blue, star-forming galaxy from the library used by 
Ilbert et al. (2008), for a total of 31 templates. 

First, we trimmed this list by discarding all templates 
that were not selected for any of the 396 sources in the 
spectroscopic sample. Next, we selected for each galaxy 
type the template which was most often used. Accord- 
ingly, the SO and the 122491 templates were singled out 
to represent passive galaxies and starburst/ULIRGs, re- 
spectively. This process reduced the original number of 
SEDs from 31 to 12. 

For a large number of objects, significant discrepancies 
between the photometry and templates occurred long- 
ward of 1 /zm. This suggests that the SEDs are not 
properly representing the mid-infrared emission compo- 
nent. The SED templates are from spectra of nearby 
galaxies where the nuclear flux is well sampled but the 



http:/ /www. iasf-milano.inaf.it/ 
letta/templates / swire_templates.html 
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host galaxy is not fully included in slits or fibers. In 
contrast, for all XMM sources at intermediate/high red- 
shift, their host galaxy sizes are comparable to the see- 
ing achievable with ground-based optical telescopes and 
our photometric data should include a larger fr action of 
the host galaxy emission (Szokolv et al.l l2"004f l. Thus, 
for the XMM sources analyzed here, the SED templates 
should very likely have a larger contribution from the 
host galaxy relative to the nuclear emission component. 

To better account for the host galaxy component, for 
each pair of AGN (type 1, type 2, QSOl, QS02) and 
host galaxy type (elliptical, spiral, starburst) hybrid tem- 
plates were constructed. Here, a pair of SEDs was nor- 
malized by the integrated fluxes in different pass bands. 
For each normalization, nine hybrids with varying ra- 
tio (90:10, 80:20, 20:80, 10:90) between both com- 
ponents were built. After re-calculating the photomet- 
ric redshifts using these hybrid SEDs, two sets of tem- 
plates provided the smallest crAz/{i+Zsp^a) number of 
outliers, rj ( defined as percentage of objects for which 
\zphot - Zspecl > 0.15(1-1- Zspec). Thcsc scts of hybrids 
are: 

• starburst/ULIRG 122491 plus TQSOl templates, both 
normahzed between 5000A and 5200A (models 19-27 in 
Tableland in Figure O; • SO plus QS02 templates, 
both normalized between 0.9/im and Ifim (models 9-17 
in Tableland in Figure O. 




0.1 1 10 100 0.1 1 10 

X(^m) X(jum) 



Fig. 3. — Visualization of the final choice of templates used in 
this paper. Numbers correspond to the IDs in Table |2] 

Lastly, a number of sources showed discrepancies be- 
tween the observed and predicted flux in the UV. This 
is likely the result of intrinsic absorption in the empiri- 
cal templates. Thus, we extended the QSO dominated 
templates 25-30 (templates with prefix "pi") into the 
UV using a. power law comp onent with spectral index 
a = -0.56 (|Scott et al.ll2004f ). Our final hbrary includes 
30 SEDs. 

4.3. Intrinsic Extinction 



TABLE 2 

Final template library used in this work. 



ID 


Namc^ 


Type 


01 


BC03 


Blue starforming"^ 


02 


SO 


SO-* 


03 


Sb 


Spiral b^ 


04 


Spi4 


Spiral c"^ 


05 


M82 


Starburst 


06 


122491 


Starburst/ULIRG^ 


07 


oey lo 


ocyioiL i.o 


08 




Spvfort 2^ 


09 


S0-90.QSO2-10 


Hybrid'* 


10 


S0-80_QSO2-20 




11 


S0-70_QSO2-30 




12 


S0-60_QSO2-40 




13 


S0-50_QSO2-50 




14 


So-40_QSO2-60 




15 


S0-30_QSO2-70 




16 


S0-20.QSO2-80 




17 


S0-10_QSO2-90 




1 s 


iVilKZOi 


oeyierL i, iD/\ijVi^ow 


1 Q 


T994Q1 QO TO'^m 10 


Hybrid'* 


90 

zu 


T99/1Q1 SO XOCJOI 90 




21 


199401-70 T0S01-3n 


„ 


22 


I22491-60_TQSOl-40 




23 


I22491-50.TQSO1-50 




24 


I22491-40.TQSO1-60 




25 


plJ22491-30.TQSOl-70 




26 


plJ22491-20.TQSOl-80 




27 


plJ22491-10.TQSOl-90 




28 


pl.QSOH 


QSO high lum.5 


29 


pLQSO 


QSO low lum.6 


30 


pLTQSOl 


QSO high IR lum.3 



To allow for the intrinsic reddening, two extinction 
laws were tested: the SMC extinction derived by Pre- 
vot et al. (1984) and the starburst galaxy extinction 
law obtained by Calzetti et al. (2000). We explored 
E(B-V) values from to 0.5 in steps of 0.05. The 
two laws provide a similar photometric redshift accuracy 
(o'Az/(i+2,pec) = 0.014), but with the SMC extinction 
law, the number of outliers is reduced by 1.1%. In prin- 
ciple, highly obscured AGN might have Ay > 1. How- 
ever, none of the fits required E(B-V) > 0.4 and 66% 
of the sources were fit with E(B-V)< 0.05. Less than 
20 sources required E(B-V)>0.35. This result is perhaps 
not totally surprising, since many of our templates are 
empirical and thus already include some extinction. 

4.4. Luminosity Prior 

Applying a luminosity prior during the SED fitting can 
be important to avoid unlikely combinations of source 
type, luminosity and redshift which can lead to wrong 
solutions for the photometric redshift. Quasars are usu- 
ally defined as point-like sources which have an absolute 
magnitude of Ms < —23. For lower luminosity AGN, the 
host galaxy contributes to the observed flux and the ob- 
ject may be classified as extended in good seeing data or, 
even more, in HST data. In the spectroscopic sample we 
did not find any point-like source with an absolute mag- 
nitude Mb fainter than —20 and thus we adopted this 
value as the lower limit on the luminosity prior. While it 
is known that low luminosity AGN may have Mb fainter 
than our limit, to use a lower threshold of Mb— -19 (as 
used by PoUetta et al. 2007 and Rowan-Robinson et al. 
2008), would increase the degeneracy between galaxy and 
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AGN templates and thus the outhers from 6.3% to 8.0%. 

Type 1 AGN and QSO can be identified first by their 
point-hke appearance in the image or by the amount of 
variabihty as defined in § [3] For the sources detected in 
the ACS/FW814 filter and brighter than FW814ab=24 
the peak surface brightness is measure d and used to sep- 
arate extended and point-like sources (jLeauthaud et al.l 
|2007() . For the 500 sources missing this parameter the 
star/galaxy separation computed on the best seeing Sub- 
aru R-band image has been used to separate extended 
and point-like sources. The sub-sample of sources with 
an extended optical counterpart includes 989 sources, 
while the sub-sample of point like sources includes 553 
members. 

We also applied the prior in luminosity to sources with 
high variability (§ [3]). For extended sources, the vari- 
ability parameter VAR can additionally be used to dis- 
entangle between host and AGN dominated sources. As 
expected, the majority of the extended sources (Figure[2l 
red distribution) exhibit only low levels of variability 
{VAR <0.25). However, a tail of higher VAR is present 
for low redshift type 1 AGN, where apparently both the 
host galaxy and the variability of the nucleus are de- 
tected. In order to identify the best threshold, we run 
the photometric redshift code using a wide range of VAR 
parameters above which the luminosity prior was applied 
also to extended objects. The prior is always used for 
point-like sources. The smallest number of outliers was 
achieved using a separation at VAR = 0.25. Accordingly, 
we found 252 extended sources in the XMM-COSMOS 
sample with VAR > 0.25. 

In the following we will refer to the point-like and to the 
variable sources (805) as "qsov" sample. To this sample 
we applied the luminosity prior. All the remaining ex- 
tended sources (737) with VAR < 0.25 form the "extnv" 
sample. No prior is applied to the latter sample. 

5. RESULTS 
5.1. Photometric redshift accuracy 

Following the procedure described in the previous sec- 
tions, we obtained photometric redshifts for all 1542 opti- 
cal counterparts of the XMM-COSMOS sample described 
in § [21 For the 396 bright spectroscopically confirmed 
galaxies we obtained a final accuracy of crAz/(i+Zsp^^} — 
0.014 and 77=4.0% outliers (see Figure [5T|) . Considering 
the "qsov" and "extnv" samples separately, the accu- 
racies are crAz/(n-2,p„)=0-012 and crA2/(n-2,p„)=0-019 
while the outliers are ri—6.3% and 7^=2.3%, respectively 
(see the first row of Table |4]) . 

Figure 15.11 shows the distribution of (zspec ~ 
Zphot) / i^+^spec) for the "extnv" and "qsov" samples. 
The histogram for the "extnv" sample is broader, not 
perfectly centered at and with a small asymmetric tail 
toward positive Az (i.e. Zphot < Zgpec)- We found that 
the scatter could be reduced for the "extnv" sample by 
adding additional normal galaxy SED templates, but this 
also increases the degeneracy and the number of outliers. 

We were also able to characterize the nature of the 
small number of outliers and to understand the reasons 
why the redshifts were not correctly assigned. First, 6 
out of the 16 outliers (37.5%) are blended with a nearby 
object, closer than l."5, while the spectroscopic follow 
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Fig. 4. — Comparison of photometric and spectroscopic red- 
shifts for the 396 sources with i\g <22.5. Empty circles represent 
sources for which the second peak in the redshift probability distri- 
bution agrees with the spectroscopic redshift. The hnes correspond 

to 1) Zphot = Zgpec (thick sohd), 2) Zphot = Zspec±0.05{l+Zspec) 

(solid) and 3) Zphot=Zspec ± 0.15(1 + Zspec) (dotted). 
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Fig. 5. — Histogram of Az/(l+z) for the "extnv" (red) and 
"qsov" (black) samples (see text for details). 



up was centered on the right counterpart, identified on 
our ACS images. In the entire XMM-COSMOS catalog 
we fiagged 34 cases as potentially problematic, due to the 
presence of a close source. Additionally, 2 of the remain- 
ing outliers miss photometry in the IRAC band (i.e. the 
source is either partially or completely in a masked area) 
and the remaining optical bands were not sufficient to 
find the correct redshift solutions. Lastly, we found that 
for strongly varying sources {VAR > 0.25) the available 
u* fluxes were likely the problem - the adopted u* mag- 
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nitudes are an average between 4 runs of observations 
taken yearly between 2004 and 2007. Adding 0.2 mag in 
quadrature to the measurement uncertainty associated 
with the u* band photometry for the outhers, 3 of the 
remaining 8 outHers would now have the correct redshift 
due to the reduced weighting of the u* band; an addi- 
tional one would have the second comparable but lower 
peak of the redshift probability distribution at the right 
redshift. We decided against systematically reducing the 
weighting of the u* band for the entire sample, as this 
would degrade the quality of the final result. 

5.2. Redshift and Templates Distribution 

Figure [5] shows the photometric redshift distribution 
for the complete XMM sample. As expected, the large 
majority of the sources at z < 1 are extended and not 
strongly time variable, while most of the point-like and 
varying sources are at z > 1. The redshift distribution 
shows significant peaks at z^^O.SS, ^^0.8, corresponding to 
spectroscopicall y confirmed galaxy oyer-densities in the 
COS MOS field fScoville erS]l2007ayFinoguenov"^^ 
[20071: iLillv eiTall 120071 : iGilfi et al.l[2008h . This suggests 
that the clustering of the normal galaxies and the AGN 
follow each other. 




Fig. 6. — Photometric redshift distribution for the entire sample 
(black line), "extnv" (red) and "qsov" (blue) samples. 



The distribution of the selected SED templates as a 
function of redshift is shown in Figure [71 Templates of 
normal galaxies (templates 1 to 6) or type 2 low luminos- 
ity AGN (e.g. Seyfert 1.8, template 7) are most frequent 
at z<2. At higher redshift, mainly QSO templates (tem- 
plates 28 to 30) have been selected. For some of the 
lower redshift galaxies, the detected X-ray emission may, 
in fact, arise from star forming regions rather than nu- 
clear activity, while at high redshift, the X-ray emission 
is predominantly from AGN. 

In Section 15.41 we show that the best fit SED tem- 
plates are generally consistent with other properties of 
the host galaxies. We note that Figure [7| clearly reveals 
that templates 14 to 17, with an increasing contribu- 



tion of the QS02 SED, were rarely selected. However, 
running the photometric redshift fitting without these 
templates resulted in a significant increase in outliers, 
indicating that for a few sources these templates are es- 
sential. In summary, we classify 464 sources as best fit by 
a non-active galaxy (templates 1-6) while 1032 are best 
fit with a template which include an AGN component 
(templates 7-30). 




10 15 20 

Templates index 



Fig. 7. — Density map of SED templates distribution as a 
function of photometric redshift. 



5.3. Optically undetected sample 

As discussed in g2 10 XMM-COSMOS sources are 
not detected in our optical broad or intermediate fil- 
ters. However, unique counterparts for these sources 
have been identified in the K-band and in the IRAC 
catalogs. These are potentially very high redshift 
sources (z> 3; Koekemoer et al. 2007, Brusa et al. 



2008), as sugg ested for similar MIPS 24 select ed objects 
(iHouck et alJ 120051: [Weedman et all l2006bl : lYan et all 
120071 : iSaiina et all l2007[ )~ However doubts on this in- 
terpretation has recently been presented by [Brand et al.l 
^00§) who observed 16 X -ray loud, optically faint 
sources with Spitzer/IRS and measured redshifts Z'^2 for 
all targets. 

Our redshift estimates for 6 of the 10 sources in 
the XMM-COSMOS optically faint sample is shown in 
Fig. 15.31 We do not report the results for the remaining 
4 sources, as strong contamination from a nearby object 
does not allow accurate photometry. In 4 of the 6 cases, 
the formal best fit redshift is higher than 4, but the red- 
shift probability distribution function (see insert in each 
figure) indicates that there are not enough constraints 
to reject a solution at lower redshift. In particular, z ~ 
2.0-2.5 would be acceptable for all of them. 

5.4. Reliability of the best fit SED templates 

It is reasonable to ask if the selected SED templates 
are consistent with other properties of the sources as re- 
vealed in the COSMOS multi-wavelength data set. For 
example, 1) are the highly varying sources mostly fit with 
a template including a type 1 AGN contribution? and 
2) were the sources likely to be "host dominated" on the 
basis of their low X-ray luminosity better fit by templates 
with little overall AGN contribution ? 

Fig. [5| clearly shows that most (^ 62%) of the strongly 
varying sources are best fit with type 1 AGN or QSO 
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l)Zph,>i= S.09 templale=28 
-2)Zph„= 7-62 template=29 



l)2phoi= 2.11 templale=l 
|2)Zph<.i" '-38 template=2 




Fig. 8. — SED fitting and piiotometric redshift estimation for the isolated 6 K band detected sources for which only 5 (or less) 
photometric points are available. The redshift probability distribution is shown in the inserts. Although the 1-cr errors (black 
horizontal bars in the insert) are so large that a solution at z~2 can not be excluded, the best formal fit redshift is greater than 
4. 
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Fig. 9. — Distribution of the best fit SED templates for all 
variable sources (black), for the point-like (blue) and the extended 
(red) sub-samples. The vertical lines indicate the separation be- 
tween templates of non-active galaxies (templates 1 to 6), type 2 
AGN or with a QS02 component (templates 7 to 17) and type 1 
AGN, QSO or with a QSO component (templates 18 to 30). The 
bulk of the sources is best fit with a template of the third group. 

dominated hybrid SEDs while ~ 25% are best fit with 
one of the non-active galaxies. Similar reassurance is 
found in the correlation of type 1, type 2 and nor- 
mal galaxies with X-ray properties like luminosity and 
hardness ratio. For example, Figure \W\ shows the 
distribution of variable (right panel) and non-variable 
(left panel) sources in the XMM-color plane. Type 
1 AGN are expected to be in a narrow region with 
HRl'^ —0.5, while typ e 2 AGN a nd normal galaxi es have 
highe r HRl values fHasin ger et al. 2001; Mainieri et al.l 
20021: iDe lla Ccca ct al.l 120041 : iHasinger et al.l 120071 : 



Cappelluti et al., ,2008. ). Our classification is well con- 
sistent with the distribution in the plane of the XMM 
sources (see caption of the figure). Note that for some 
sources, the X-ray flux is obtained by stacking up to 9 
pointings taken over 3 years. Thus, the determination of 
the HR values can be effected by spectral variability. 

The 2-10 keV X-ray luminosity is commonly used to 
distinguish AGN (Lx > 10^^ erg/s) from starburst or 
normal galaxies (Lx < lO'*^ erg/s). As we show in Fig- 
ure llli the large majority of low-luminosity sources are 
non- variable and a host dominated template was selected 
for most of them. A large fraction of high-luminosity 
objects is instead varying and the best fit is obtained 
with type 1 and QSO SEDs. A number of sources with 
(Lx > lO'^^ erg/s) have VAR < 0.25 and are fit by a 
normal galaxy template. This is consistent with cases in 
which the optical nuclear activity is completely hidden 
by the host galaxy. 

5.5. Reliability of the photometric redshift for faint 
sources 

Recently, the spectroscopic redshift for 46 faint sources 
22.5 <i^5 <24.5 (average value i^^ ^23.3) became 
available. They are part of a Keck II spectroscopic 
follow-up of 24/im selected sources (Kartaltepe et al. 



2008) and part of zCOSMOS-faint (Lilly et al. 2008). 
Fig. 15.5! (left panel) shows the comparison between the 
photometric and spectroscopic redshifts for this sample. 
As these sources are faint, the photometric uncertain- 
ties are larger and we should expect lower photometric 
redshift accuracy. In fact, the accuracy become worse 
(although still comparable to the typical accuracy of pho- 
tometric redshift for non-active galaxies) by a factor of 
1.5, from 0.014 to 0.023, while the number of outliers 
increases from 4.0% (for the bright sample) to 10.9%. 
However, two of the five outliers (open circles) do have 
a secondary maximum in the redshift probability distri- 
bution at the spectroscopic value. 

As the reliability of the photometric redshift based on 
the comparison with a spectroscopic sample is limited 
to the redshift and magnitude ranges of the the sample, 
we used the la uncertainty (as derived from the redshift 
probability distribution function) to test the accuracy 
of the photometric redshift for all the AMM-COSMOS 
sources. Figure l575l (right panel) shows the la uncertain- 
ties in the photometric redshift as a function of magni- 
tude for sources with 0.2<z<1.25 and 1.25<z< 5. Unlike 
normal galaxies (see, for example, Figure 9 in Ilbert et al. 
2008), the uncertainties for the two redshift ranges be- 
have in the same way and up to ij^^ ^24 remain almost 
constant, of the order of 0.02. Only beyond ij^^ >24 
the uncertainties increase and reach an amplitude of 0.2 



at 1 



AB 



^26. This is expected as AGN show prominent 
emission lines, all over the spectrum. In particular, be- 
side the well known optical lines (Hbeta, [OIII], Haipha, 
etc.) two strong emission lines in the blue part of the 
spectrum (CIV[1549A] and MgII[2798A]) are available 
at high redshift. The COSMOS survey, dense with in- 
termediate band photometry, allows us to identify these 
lines and thus provides a reliable photometric redshift 
estimate over the full redshift range. The reliability of 
our photometric redshift is diminished only by the larger 
photometric errors at faint magnitudes. 

5.6. Reliability of star/galaxy separation 

In addition to galaxy and AGN templates, stel- 
lar templates have been used to identify stars in the 
AMM-COSMOS sampl e We constructed a s pectral li- 
brary of low mass stars (|Chabrier et al. I [20001 ) and sub- 
dwarf O and B sta rs, white dwarfs and binary systems 
(jBixler et al.l I1991D . The reliability of the star/galaxy 
separation was tested on the sub-sample with spectro- 
scopic identification. We defined as star each point-like 
source with 

(2) 



2 

Xstar 



and the remaining as extra-galactic. In this way we 
classify 46 sources as stars, 28 of which are already 
spectroscopically confirmed. In summary, 439 (out of 
396-1-46=442) galaxies and 24 (out of 28) stars were clas- 
sified correctly. The completeness in identifying galaxies 
is 439/442= 0.99 and 24/28= 0.86 for the stars. Two of 
the three misclassified galaxies were missing IRAC pho- 
tometry and thus the possibility to disentangle between 
star and galaxies was limited. 

6. RELEASE OF THE PHOTOMETRIC REDSHIFT 
CATALOG 



Photometric redshifts for AGN 



11 



0.5 



S 



-0.5 



1 1 

- 


1 1 1 1 1 1 1 1 1 1 


1 

- 








- 




; 




- 








- N„=; 

N =C 













0.5 



^ 



-0.5 



3 2 10 




-1 



-0.5 



0.5 



-1 



-0.5 



0.5 



HR2 



HR2 



Fig. 10.— Left: X-ray hardness ratios of all XMM-COSMOS non varying sources (V AR < 0.25) and with AHRK 0.25 and 



AHR2< 0.25. The hardness ratios are defined as HR1=(F( 



2-10keV)-i^ (0.5-2kcV) 



)/(F( 



2-lOkeV) 



+F(o.5-2keV)) and HR2=(F( 



5-lOkcV)- 




lg(l,)< 42 

■lg(lT,)< 42; VAR<0.25 



F(2-i0kGV) )/(F(5— iOkeV)+F(2-iokcV) )• The best fit templates are color coded as follows: blue filled circles (Type 1 and QSOl hybrids — 
templates #18-30), red open triangles (Type 2 and QS02 hybrids - templates #7—17) and green open circles (normal galaxies — templates 
#1—6). The solid lines represent lines of constant Njj and the dotted lines are lines of equal photon index. Right: same as left, for sources 
with VAR > 0.25. 

while the number of outhers can be decreased with the 
appropriate SED template library and the correction for 
variability, a good photometric coverage is essential in 
order to improve the accuracy of the photometric red- 
shifts. 



7.1. Relevance of SED templates library 

An X-ray selected sample of sources consists either of 
low redshift star-forming galaxies or AGN and QSOs and 
the choice of templates must represent these populations. 
To verify this, we first compared our photometric red- 
shifts with those obtained for the same sources by Ilbert 
et al. (2008) using the same photometric data set but 
only normal galaxy SED templates. The result is shown 
in the first row of Table [31 For the "etxnv" sample, 
where the sources are galaxy dominated, the library of 
Ilbert et al. (2008) provides an excellent accuracy al- 
though the number of outliers is large. For the "qsov" 
sample dominated by AGN and QSOs, the lack of appro- 
priate templates leads to an expected failure^^. 

The addition of AGN and QSO templates in the li- 
brary does not automatically imply an improvement in 
the photometric redshifts. We verified this by comput- 
ing the photometric redshift using the library of Polletta 
et al. (2007) which includes normal galaxies, AGN and 
QSOs. The result of the test is shown in the second 
row of Table [31 In this case, the accuracy is very good 
and comparable with the one obtained with our library. 
However, the number of outliers remain high. We then 
conclude that our library, constructed as described in 
§ I4.2i can be considered as representative of a typical 
population of X-ray emitting galaxies. 
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Fig. 11. — SED template distribution for all sources split on 
the basis of their X-ray luminosity (bottom: hx < 10^^ erg/s, top: 
Lx > 10^^ erg/s) and variability (solid line: all, dashed: VAR < 
0.25). Practically all sources with Ljf < 10*'^ erg/s are non- variable 
(bottom panel). 

An electronic version of the complete photometric red- 
shift and classification catalog of the XMM-COSMOS 
sample is available. For each source the coordinates, the 
photometric redshift and the classification in star/galaxy 
is provided. A more comprehensive catalog is available 
by request'^". 

7. DISCUSSION 

In this section we analyze the impact of each contri- 
bution (new template library, correction for variability, 
photometric coverage) on our final result. We show that 

ms@astro.caltech.edu 



31 In Ilbert et al. (2008) the new publicly available COSMOS 
photometric redshift catalog is presented. This catalog supersedes 
all previous ones (Mobasher et al. 2007, Ilbert et al. 2006) and is 
composed of 600,000 sources with ij^^ <26 and not detected in X- 
rays. In order to decrease the risk of degeneracies for this sample of 
normal galaxies, no AGN/QSOs templates were used to compute 
the photometric redshifts. 
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spectro-z i'^B 

Fig. 12. — LeftiSame as Figure 4 for the 46 sources with ij[^ >22.5 with spectroscopic redshifts, which has been used for an a posteriori 
test. Right; 1-cr error for the Zp/^at estimate as a function of the apparent magnitude in the redshift range 0.2 < Zpf^^t ^ 1-25 (blue) and 
1.25 < Zphot < 5 (red). 



7.2. Relevance of correction for variability 

In the previous section we have shown the importance 
of the hbrary to increasing the accuracy and to reducing 
the number of outhers. The comparison has been per- 
formed on the photometric data set corrected for vari- 
abihty; thus, it is natural to ask how much the improved 
performances shown by our hbrary are depending on it. 
In the third row of Table [3l our library has been used 
on the photometric data set without any correction for 
variability. While the accuracy of the photometric red- 
shift did not change, the number of outliers is increased 
by a factor of 2 in both "extnv" and "qsov" sub-samples, 
demonstrating the importance of the correction for vari- 
ability. 

7.3. Relevance of spectral coverage 

To quantify the benefits of spectral coverage, we split 
the available photometric data into 5 sub-sets {GALEX, 
optical broad band, optical intermediate band, near- 
infrared, and IRAC). Next, we derived the photomet- 
ric redshifts for the "extnv" and "qsov" samples using 
different combinations of these sub-sets and compared 
the results in terms of accuracy and number of outliers. 
We find that for all the combinations the accuracy is 
o'Az/d+z, ec) < 0.07, which is good enough for most AGN 
studies. X summary is given in Tabled 

For the "extnv" sample, the by far best result was 
achieved using the entire spectral coverage. This is not 
surprising and is also discussed for normal galaxies in 
COSMOS (Ilbert et al. 2008). The 30 band photometry 
of COSMOS essentially resembles a low-resolution spec- 
trum and a representative SED template library can be 
fit with high reliability. 

Table 13] shows also the importance of the intermediate 
optical filters which allow the precise localization of spec- 
tral features (e.g., emission and absorption lines, Balmer 
break). The elimination of these filters produces a sig- 
nificantly worse accuracy (by a factor of 3 for both the 
"qsov" and the "extnv" samples). At the same time the 
number of outliers increases by a factor of 3 for both 



samples. 

In contrast, the removal of the broad optical bands 
photometry when the intermediate bands are available, 
does not change significantly the accuracy of the photo- 
metric redshift, neither for the "qsov" nor for the "extnv" 
samples. However, in both cases the number of outliers 
increases dramatically. The same is happening when the 
IRAC photometry is not consider ed for the " extnv" sam- 
ple. A s already pointed out by iRowan-Robinson et al.l 
(|200l the presence of IRAC photometry is crucial in 
suppressing the number of outliers. 

8. CONCLUSIONS 

In this paper we presented the photometric redshifts 
and classifications for 1542 AMM-COSMOS sources with 
optical counterpart. Using SED matching to up to 30 
photometric bands we identified 46 Galactic foreground 
stars, 464 sources best fit with a non-active galaxies 
template and 1032 best fit by a template which has 
an AGN contribution. A high reliability in recogniz- 
ing stars from galaxies and AGN was demonstrated us- 
ing a spectroscopically identified sub-set of 470 sources. 
Overall, we reach a photometric redshift accuracy of 
o'Az/(i+z,peo) = 0-014 (i^B < 22.5) with 4.0% outliers. 
However, the most important result is perhaps the ac- 
curacy for the sub-sample of AGN dominated sources 
('^Az/(i+2bpcc) ~ 0.011) which marks an significant im- 
provement in photometric redshift estimation for this 
kind of sources. The number of outliers, however, is still 
non-negligible (6.3% when considering the "qsov" sample 
alone). 

The excellent photometric redshifts and classifications 
have been made possible by a number of factors. Here, 
the most important was likely the large spectral cover- 
age from UV to 8 ^m, densely sampled with broad, in- 
termediate and narrow filters. The 25 optical bands es- 
sentially mimic low resolution spectra and are powerful 
in constraining the positions of spectral features. The 
near/mid-infrared coverage provides the spectral range 
required to better distinguish stars from galaxies and 
helps to constrain SED templates, especially for passive 
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TABLE 3 

Relevance of variability correction and template library 





ex 

O'A^/Cl+z) 


tnv 
outliers (%) 


q 

O'A^/CH-;!) 


30V 

outliers (%) 


IlbertOS^ 


0.014 


14.0 


0.419 


71.4 


PollettaOy^ + var. corr. 


0.020 


12.8 


0.014 


15.0 


SalvatoOS"^, no var. corr. 


0.017 


4.1 


0.013 


11.6 


SalvatoOS^ + var. corr. 


0.019 


2.3 


0.012 


6.3 



TABLE 4 

Relevance of spectral coverage and photometric resolution for the variability corrected bright spectroscopic sample. 



Sample 


GALEX 


Optical 
Broad Band 


Optical 
IB & NB 


NIR 


IRAC 


Reliat 
outliers (%) 


)ility 

0"A^/(1 + b) 


qsov 
extnv 


^/ 










6.3 
2.3 


0.012 
0.019 


qsov 
extnv 








^/ 


^/ 


19.4 
6.8 


0.063 
0.056 


qsov 
extnv 


V 










6.9 
5.9 


0.011 
0.019 


qsov 
extnv 












14.3 
6.8 


0.012 
0.021 



galaxies. 

Also important was the recognition of and correction 

for variability. This is especially relevant for AGN dom- 
inated sources, for which photometry collected over a 
large time span may dramatically vary from a snapshot 
SED. We applied a correction for variability by re-scaling 
the optical photometry to a common epoch close to that 
of IRAC observations. 

The correction decreased the number of outliers by a 
factor of 2. The significant improvement obtained in our 
sample with this correction demonstrates the importance 
of quasi-simultaneous observations for projects aiming 
for good photometric redshift estimates of AGN. Future 
surveys will have to consider this point very carefully. 
Deep imaging in a common filter in each observing epoch 
would provide the possibility of a valuable first order 
correction. 

Finally, we used a revised set of SED templates to fully 
represent the entire range of different contributions of 
active galactic nuclei and their host galaxies. Empirical 



SED templates are often obtained from objects in the 
nearby Universe where the nuclear and extended fluxes 
arc easily resolved. Their application to high-redshift 
objects, where the two contributions are mixed together, 
may not always be appropriate. Our hybrid templates 
already showed a significant improvement, although fur- 
ther advancements are required. 

The catalog of XMM-COSMOS selected sources with 
photometric redshifts and classifications, together with 
the X-ray properties (Brusa et al. 2008) offers an un- 
precedented possibility to study the luminosity function 
and redshift distribution of AGN. In combination with 
the immense multi-band COSMOS data set and the new 
photometric redshifts for optically selected sources (II- 
bert et al. 2008), studies of the environment of AGN 
will also allow to probe AGN fuelling models, clustering 
and evolution. 
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